Macrophage mediated mesoscale brain mechanical homeostasis mechanically imaged via optical tweezers and Brillouin microscopy in vivo.

Tissues are active materials where epithelial turnover, immune surveillance, and remodeling of stromal cells such as macrophages all regulate form and function. Scattering modalities such as Brillouin microscopy (BM) can non-invasively access mechanical signatures at GHz. However, our traditional understanding of tissue material properties is derived mainly from modalities which probe mechanical properties at different frequencies. Thus, reconciling measurements amongst these modalities remains an active area. Here, we compare optical tweezer active microrheology (OT-AMR) and Brillouin microscopy (BM) to longitudinally map brain development in the larval zebrafish. We determine that each measurement is able to detect a mechanical signature linked to functional units of the brain. We demonstrate that the corrected BM-Longitudinal modulus using a density factor correlates well with OT-AMR storage modulus at lower frequencies. We also show that the brain tissue mechanical properties are dependent on both the neuronal architecture and the presence of macrophages. Moreover, the BM technique is able to delineate the contributions to mechanical properties of the macrophage from that due to colony stimulating factor 1 receptor (CSF1R) mediated stromal remodeling. Here, our data suggest that macrophage remodeling is instrumental in the maintenance of tissue mechanical homeostasis during development. Moreover, the strong agreement between the OT-AM and BM further demonstrates that scattering-based technique is sensitive to both large and minute structural modification in vivo.

Genomic landscape of patients with germline RUNX1 variants and familial platelet disorder with myeloid malignancy.

ABSTRACT: Familial platelet disorder with associated myeloid malignancies (FPDMM) is caused by germline RUNX1 mutations and characterized by thrombocytopenia and increased risk of hematologic malignancies. We recently launched a longitudinal natural history study for patients with FPDMM. Among 27 families with research genomic data by the end of 2021, 26 different germline RUNX1 variants were detected. Besides missense mutations enriched in Runt homology domain and loss-of-function mutations distributed throughout the gene, splice-region mutations and large deletions were detected in 6 and 7 families, respectively. In 25 of 51 (49%) patients without hematologic malignancy, somatic mutations were detected in at least 1 of the clonal hematopoiesis of indeterminate potential (CHIP) genes or acute myeloid leukemia (AML) driver genes. BCOR was the most frequently mutated gene (in 9 patients), and multiple BCOR mutations were identified in 4 patients. Mutations in 6 other CHIP- or AML-driver genes (TET2, DNMT3A, KRAS, LRP1B, IDH1, and KMT2C) were also found in ≥2 patients without hematologic malignancy. Moreover, 3 unrelated patients (1 with myeloid malignancy) carried somatic mutations in NFE2, which regulates erythroid and megakaryocytic differentiation. Sequential sequencing data from 19 patients demonstrated dynamic changes of somatic mutations over time, and stable clones were more frequently found in older adult patients. In summary, there are diverse types of germline RUNX1 mutations and high frequency of somatic mutations related to clonal hematopoiesis in patients with FPDMM. Monitoring changes in somatic mutations and clinical manifestations prospectively may reveal mechanisms for malignant progression and inform clinical management. This trial was registered at www.clinicaltrials.gov as #NCT03854318.

Fluorescent PCR-based Screening Methods for Precise Knock-in of Small DNA Fragments and Point Mutations in Zebrafish.

Generation of zebrafish (Danio rerio) models with targeted insertion of epitope tags and point mutations is highly desirable for functional genomics and disease modeling studies. Currently, CRISPR/Cas9-mediated knock-in is the method of choice for insertion of exogeneous sequences by providing a repair template for homology-directed repair (HDR). A major hurdle in generating knock-in models is the labor and cost involved in screening of injected fish to identify the precise knock-in events due to low efficiency of the HDR pathway in zebrafish. Thus, we developed fluorescent PCR-based high-throughput screening methods for precise knock-in of epitope tags and point mutations in zebrafish. Here, we provide a step-by-step guide that describes selection of an active sgRNA near the intended knock-in site, design of single-stranded oligonucleotide (ssODN) templates for HDR, quick validation of somatic knock-in using injected embryos, and screening for germline transmission of precise knock-in events to establish stable lines. Our screening method relies on the size-based separation of all fragments in an amplicon by fluorescent PCR and capillary electrophoresis, thus providing a robust and cost-effective strategy. Although we present the use of this protocol for insertion of epitope tags and point mutations, it can be used for insertion of any small DNA fragments (e.g., LoxP sites, in-frame codons). Furthermore, the screening strategy described here can be used to screen for precise knock-in of small DNA sequences in any model system, as PCR amplification of the target region is its only requirement. Key features This protocol expands the use of fluorescent PCR and CRISPR-STAT for screening of precise knock-in of small insertions and point mutations in zebrafish. Allows validation of selected sgRNA and HDR template within two weeks by somatic knock-in screening. Allows robust screening of point mutations by combining restriction digest with CRISPR-STAT. Graphical overview Overview of the three-phase knock-in pipeline in zebrafish (created with BioRender.com).

Somatic mutational landscape of hereditary hematopoietic malignancies caused by germline variants in RUNX1, GATA2, and DDX41.

Individuals with germ line variants associated with hereditary hematopoietic malignancies (HHMs) have a highly variable risk for leukemogenesis. Gaps in our understanding of premalignant states in HHMs have hampered efforts to design effective clinical surveillance programs, provide personalized preemptive treatments, and inform appropriate counseling for patients. We used the largest known comparative international cohort of germline RUNX1, GATA2, or DDX41 variant carriers without and with hematopoietic malignancies (HMs) to identify patterns of genetic drivers that are unique to each HHM syndrome before and after leukemogenesis. These patterns included striking heterogeneity in rates of early-onset clonal hematopoiesis (CH), with a high prevalence of CH in RUNX1 and GATA2 variant carriers who did not have malignancies (carriers-without HM). We observed a paucity of CH in DDX41 carriers-without HM. In RUNX1 carriers-without HM with CH, we detected variants in TET2, PHF6, and, most frequently, BCOR. These genes were recurrently mutated in RUNX1-driven malignancies, suggesting CH is a direct precursor to malignancy in RUNX1-driven HHMs. Leukemogenesis in RUNX1 and DDX41 carriers was often driven by second hits in RUNX1 and DDX41, respectively. This study may inform the development of HHM-specific clinical trials and gene-specific approaches to clinical monitoring. For example, trials investigating the potential benefits of monitoring DDX41 carriers-without HM for low-frequency second hits in DDX41 may now be beneficial. Similarly, trials monitoring carriers-without HM with RUNX1 germ line variants for the acquisition of somatic variants in BCOR, PHF6, and TET2 and second hits in RUNX1 are warranted.

A novel de novo TP63 mutation in whole-exome sequencing of a Syrian family with Oral cleft and ectrodactyly.

BACKGROUND: Oral clefts and ectrodactyly are common, heterogeneous birth defects. We performed whole-exome sequencing (WES) analysis in a Syrian family. The proband presented with both orofacial clefting and ectrodactyly but not ectodermal dysplasia as typically seen in ectrodactyly, ectodermal dysplasia, and cleft lip/palate syndrome-3. A paternal uncle with only an oral cleft was deceased and unavailable for analysis. METHODS: Variant annotation, Mendelian inconsistencies, and novel variants in known cleft genes were examined. Candidate variants were validated using Sanger sequencing, and pathogenicity assessed by knocking out the tp63 gene in zebrafish to evaluate its role during zebrafish development. RESULTS: Twenty-eight candidate de novo events were identified, one of which is in a known oral cleft and ectrodactyly gene, TP63 (c.956G > T, p.Arg319Leu), and confirmed by Sanger sequencing. CONCLUSION: TP63 mutations are associated with multiple autosomal dominant orofacial clefting and limb malformation disorders. The p.Arg319Leu mutation seen in this patient is de novo but also novel. Two known mutations in the same codon (c.956G > A, p.(Arg319His; rs121908839, c.955C > T), p.Arg319Cys) cause ectrodactyly, providing evidence that mutating this codon is deleterious. While this TP63 mutation is the best candidate for the patient's clinical presentation, whether it is responsible for the entire phenotype is unclear. Generation and characterization of tp63 knockout zebrafish showed necrosis and rupture of the head at 3 days post-fertilization (dpf). The embryonic phenotype could not be rescued by injection of zebrafish or human messenger RNA (mRNA). Further functional analysis is needed to determine what proportion of the phenotype is due to this mutation.

Genomic Landscape of Patients with Germline RUNX1 Variants and Familial Platelet Disorder with Myeloid Malignancy.

Germline RUNX1 mutations lead to familial platelet disorder with associated myeloid malignancies (FPDMM), which is characterized by thrombocytopenia and a life-long risk (35-45%) of hematological malignancies. We recently launched a longitudinal natural history study for patients with FPDMM at the NIH Clinical Center. Among 29 families with research genomic data, 28 different germline RUNX1 variants were detected. Besides missense mutations enriched in Runt homology domain and loss-of-function mutations distributed throughout the gene, splice-region mutations and large deletions were detected in 6 and 7 families, respectively. In 24 of 54 (44.4%) non-malignant patients, somatic mutations were detected in at least one of the clonal hematopoiesis of indeterminate potential (CHIP) genes or acute myeloid leukemia (AML) driver genes. BCOR was the most frequently mutated gene (in 9 patients), and multiple BCOR mutations were identified in 4 patients. Mutations in 7 other CHIP or AML driver genes ( DNMT3A, TET2, NRAS, SETBP1, SF3B1, KMT2C , and LRP1B ) were also found in more than one non-malignant patient. Moreover, three unrelated patients (one with myeloid malignancy) carried somatic mutations in NFE2 , which regulates erythroid and megakaryocytic differentiation. Sequential sequencing data from 19 patients demonstrated dynamic changes of somatic mutations over time, and stable clones were more frequently found in elderly patients. In summary, there are diverse types of germline RUNX1 mutations and high frequency of somatic mutations related to clonal hematopoiesis in patients with FPDMM. Monitoring dynamic changes of somatic mutations prospectively will benefit patients’ clinical management and reveal mechanisms for progression to myeloid malignancies. Key Points: Comprehensive genomic profile of patients with FPDMM with germline RUNX1 mutations. Rising clonal hematopoiesis related secondary mutations that may lead to myeloid malignancies.

Organ specific microenvironmental MR1 expression in cutaneous melanoma.

The microenvironment is an important regulator of intertumoral trafficking and activity of immune cells. Understanding how the immune system can be tailored to maintain anti-tumor killing responses in metastatic disease remains an important goal. Thus, immune mediated eradication of metastasis requires the consideration of organ specific microenvironmental cues. Using a xenograft model of melanoma metastasis in adult zebrafish, we perturbed the dynamic balance between the infiltrating immune cells in the metastatic setting using a suite of different transgenic zebrafish. We employed intravital imaging coupled with metabolism imaging (FLIM) to visualize and map the organ specific metabolism with near simultaneity in multiple metastatic lesions. Of all the MHC complexes examined for brain and skeletal metastases, we determined that there is an organ specific expression of mhc1uba (human ortholog, MR1) for both the melanoma cells and the resident and infiltrating immune cells. Specifically, immune clusters did not express mhc1uba in brain metastatic lesions in immune competent fish. Finally, the differential immune response drove organ specific metabolism where tumor glycolysis was increased in brain metastases compared to skeletal and parental lines as measured using fluorescence lifetime imaging microscopy (FLIM). As MR1 belongs to the MHC class I molecules and is a target of immunotherapeutic drugs, we believe that our data presents an opportunity to understand the relationship between organ specific tumor metabolism and drug efficacy in the metastatic setting.

A robust pipeline for efficient knock-in of point mutations and epitope tags in zebrafish using fluorescent PCR based screening.

BACKGROUND: Genome editing using CRISPR/Cas9 has become a powerful tool in zebrafish to generate targeted gene knockouts models. However, its use for targeted knock-in remains challenging due to inefficient homology directed repair (HDR) pathway in zebrafish, highlighting the need for efficient and cost-effective screening methods.  RESULTS: Here, we present our fluorescent PCR and capillary electrophoresis based screening approach for knock-in using a single-stranded oligodeoxynucleotide donor (ssODN) as a repair template for the targeted insertion of epitope tags, or single nucleotide changes to recapitulate pathogenic human alleles. For the insertion of epitope tags, we took advantage of the expected change in size of the PCR product. For point mutations, we combined fluorescent PCR with restriction fragment length polymorphism (RFLP) analysis to distinguish the fish with the knock-in allele. As a proof-of-principle, we present our data on the generation of fish lines with insertion of a FLAG tag at the tcnba locus, an HA tag at the gata2b locus, and a point mutation observed in Gaucher disease patients in the gba gene. Despite the low number of germline transmitting founders (1-5%), combining our screening methods with prioritization of founder fish by fin biopsies allowed us to establish stable knock-in lines by screening 12 or less fish per gene. CONCLUSIONS: We have established a robust pipeline for the generation of zebrafish models with precise integration of small DNA sequences and point mutations at the desired sites in the genome. Our screening method is very efficient and easy to implement as it is PCR-based and only requires access to a capillary sequencer.

Zrsr2 Is Essential for the Embryonic Development and Splicing of Minor Introns in RNA and Protein Processing Genes in Zebrafish.

ZRSR2 (zinc finger CCCH-type, RNA binding motif and serine/arginine rich 2) is an essential splicing factor involved in 3' splice-site recognition as a component of both the major and minor spliceosomes that mediate the splicing of U2-type (major) and U12-type (minor) introns, respectively. Studies of ZRSR2-depleted cell lines and ZRSR2-mutated patient samples revealed its essential role in the U12-dependent minor spliceosome. However, the role of ZRSR2 during embryonic development is not clear, as its function is compensated for by Zrsr1 in mice. Here, we utilized the zebrafish model to investigate the role of zrsr2 during embryonic development. Using CRISPR/Cas9 technology, we generated a zrsr2-knockout zebrafish line, termed zrsr2hg129/hg129 (p.Trp167Argfs*9) and examined embryo development in the homozygous mutant embryos. zrsr2hg129/hg129 embryos displayed multiple developmental defects starting at 4 days post fertilization (dpf) and died after 8 dpf, suggesting that proper Zrsr2 function is required during embryonic development. The global transcriptome analysis of 3 dpf zrsr2hg129/hg129 embryos revealed that the loss of Zrsr2 results in the downregulation of essential metabolic pathways and the aberrant retention of minor introns in about one-third of all minor intron-containing genes in zebrafish. Overall, our study has demonstrated that the role of Zrsr2 as a component of the minor spliceosome is conserved and critical for proper embryonic development in zebrafish.

Dysregulated myosin in Hermansky-Pudlak syndrome lung fibroblasts is associated with increased cell motility.

Hermansky-Pudlak syndrome (HPS) is an autosomal recessive disorder characterized by improper biogenesis of lysosome-related organelles (LROs). Lung fibrosis is the leading cause of death among adults with HPS-1 and HPS-4 genetic types, which are associated with defects in the biogenesis of lysosome-related organelles complex-3 (BLOC-3), a guanine exchange factor (GEF) for a small GTPase, Rab32. LROs are not ubiquitously present in all cell types, and specific cells utilize LROs to accomplish dedicated functions. Fibroblasts are not known to contain LROs, and the function of BLOC-3 in fibroblasts is unclear. Here, we report that lung fibroblasts isolated from patients with HPS-1 have increased migration capacity. Silencing HPS-1 in normal lung fibroblasts similarly leads to increased migration. We also show that the increased migration is driven by elevated levels of Myosin IIB. Silencing HPS1 or RAB32 in normal lung fibroblasts leads to increased MYOSIN IIB levels. MYOSIN IIB is downstream of p38-MAPK, which is a known target of angiotensin receptor signaling. Treatment with losartan, an angiotensin receptor inhibitor, decreases MYOSIN IIB levels and impedes HPS lung fibroblast migration in vitro. Furthermore, pharmacologic inhibition of angiotensin receptor with losartan seemed to decrease migration of HPS lung fibroblasts in vivo in a zebrafish xenotransplantation model. Taken together, we demonstrate that BLOC-3 plays an important role in MYOSIN IIB regulation within lung fibroblasts and contributes to fibroblast migration.

A Comprehensive Review of Indel Detection Methods for Identification of Zebrafish Knockout Mutants Generated by Genome-Editing Nucleases.

The use of zebrafish in functional genomics and disease modeling has become popular due to the ease of targeted mutagenesis with genome editing nucleases, i.e., zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs), and clustered regularly interspaced short palindromic repeats/Cas9 (CRISPR/Cas9). These nucleases, specifically CRISPR/Cas9, are routinely used to generate gene knockout mutants by causing a double stranded break at the desired site in the target gene and selecting for frameshift insertions or deletions (indels) caused by the errors during the repair process. Thus, a variety of methods have been developed to identify fish with indels during the process of mutant generation and phenotypic analysis. These methods range from PCR and gel-based low-throughput methods to high-throughput methods requiring specific reagents and/or equipment. Here, we provide a comprehensive review of currently used indel detection methods in zebrafish. By discussing the molecular basis for each method as well as their pros and cons, we hope that this review will serve as a comprehensive resource for zebrafish researchers, allowing them to choose the most appropriate method depending upon their budget, access to required equipment and the throughput needs of the projects.

Redundant mechanisms driven independently by RUNX1 and GATA2 for hematopoietic development.

RUNX1 is essential for the generation of hematopoietic stem cells (HSCs). Runx1-null mouse embryos lack definitive hematopoiesis and die in mid-gestation. However, although zebrafish embryos with a runx1 W84X mutation have defects in early definitive hematopoiesis, some runx1W84X/W84X embryos can develop to fertile adults with blood cells of multilineages, raising the possibility that HSCs can emerge without RUNX1. Here, using 3 new zebrafish runx1-/- lines, we uncovered the compensatory mechanism for runx1-independent hematopoiesis. We show that, in the absence of a functional runx1, a cd41-green fluorescent protein (GFP)+ population of hematopoietic precursors still emerge from the hemogenic endothelium and can colonize the hematopoietic tissues of the mutant embryos. Single-cell RNA sequencing of the cd41-GFP+ cells identified a set of runx1-/--specific signature genes during hematopoiesis. Significantly, gata2b, which normally acts upstream of runx1 for the generation of HSCs, was increased in the cd41-GFP+ cells in runx1-/- embryos. Interestingly, genetic inactivation of both gata2b and its paralog gata2a did not affect hematopoiesis. However, knocking out runx1 and any 3 of the 4 alleles of gata2a and gata2b abolished definitive hematopoiesis. Gata2 expression was also upregulated in hematopoietic cells in Runx1-/- mice, suggesting the compensatory mechanism is conserved. Our findings indicate that RUNX1 and GATA2 serve redundant roles for HSC production, acting as each other's safeguard.

Large-scale generation and phenotypic characterization of zebrafish CRISPR mutants of DNA repair genes.

A systematic knowledge of the roles of DNA repair genes at the level of the organism has been limited due to the lack of appropriate experimental approaches using animal model systems. Zebrafish has become a powerful vertebrate genetic model system with availability due to the ease of genome editing and large-scale phenotype screening. Here, we generated zebrafish mutants for 32 DNA repair and replication genes through multiplexed CRISPR/Cas9-mediated mutagenesis. Large-scale phenotypic characterization of our mutant collection revealed that three genes (atad5a, ddb1, pcna) are essential for proper embryonic development and hematopoiesis; seven genes (apex1, atrip, ino80, mre11a, shfm1, telo2, wrn) are required for growth and development during juvenile stage and six genes (blm, brca2, fanci, rad51, rad54l, rtel1) play critical roles in sex development. Furthermore, mutation in six genes (atad5a, brca2, polk, rad51, shfm1, xrcc1) displayed hypersensitivity to DNA damage agents. Our zebrafish mutant collection provides a unique resource for understanding of the roles of DNA repair genes at the organismal level.

Compound heterozygous KCTD7 variants in progressive myoclonus epilepsy.

KCTD7 is a member of the potassium channel tetramerization domain-containing protein family and has been associated with progressive myoclonic epilepsy (PME), characterized by myoclonus, epilepsy, and neurological deterioration. Here we report four affected individuals from two unrelated families in which we identified KCTD7 compound heterozygous single nucleotide variants through exome sequencing. RNAseq was used to detect a non-annotated splicing junction created by a synonymous variant in the second family. Whole-cell patch-clamp analysis of neuroblastoma cells overexpressing the patients' variant alleles demonstrated aberrant potassium regulation. While all four patients experienced many of the common clinical features of PME, they also showed variable phenotypes not previously reported, including dysautonomia, brain pathology findings including a significantly reduced thalamus, and the lack of myoclonic seizures. To gain further insight into the pathogenesis of the disorder, zinc finger nucleases were used to generate kctd7 knockout zebrafish. Kctd7 homozygous mutants showed global dysregulation of gene expression and increased transcription of c-fos, which has previously been correlated with seizure activity in animal models. Together these findings expand the known phenotypic spectrum of KCTD7-associated PME, report a new animal model for future studies, and contribute valuable insights into the disease.

Clinical and genomic analysis of a large Chinese family with familial cortical myoclonic tremor with epilepsy and SAMD12 intronic repeat expansion.

Objective: Our goal was to perform detailed clinical and genomic analysis of a large multigenerational Chinese family with 21 individuals showing symptoms of Familial Cortical Myoclonic Tremor with Epilepsy (FCMTE) that we have followed for over 20 years. Methods: Patients were subjected to clinical evaluation, routine EEG, and structural magnetic resonance imaging. Whole exome sequencing, repeat-primed PCR, long-range PCR, and PacBio sequencing were performed to characterize the disease-causing mutation in this family. Results: All evaluated patients manifested adult-onset seizures and presented with progressive myoclonic postural tremors starting after the third or fourth decade of life. Seizures typically diminished markedly in frequency with implementation of antiseizure medications but did not completely cease. The electroencephalogram of affected individuals showed generalized or multifocal spikes and slow wave complexes. An expansion of TTTTA motifs with addition of TTTCA motifs in intron 4 of SAMD12 was identified to segregate with the disease phenotype in this family. Furthermore, we found that the mutant allele is unstable and can undergo both contraction and expansion by changes in the number of repeat motifs each time it is passed to the next generation. The size of mutant allele varied from 5 to 5.5 kb with 549-603 copies of TTTTA and 287-343 copies of TTTCA repeat motifs in this family. Significance: Our study provides a detailed description of clinical progression of FCMTE symptoms and its management with antiseizure medications. Our method of repeat analysis by PacBio sequencing of long-range PCR products does not require high-quality DNA and hence can be easily applied to other families to elucidate any correlation between the repeat size and phenotypic variables, such as, age of onset, and severity of symptoms.

Functional analysis of Sonic Hedgehog variants associated with holoprosencephaly in humans using a CRISPR/Cas9 zebrafish model.

Genetic variation in the highly conserved Sonic Hedgehog (SHH) gene is one of the most common genetic causes for the malformations of the brain and face in humans described as the holoprosencephaly clinical spectrum. However, only a minor fraction of known SHH variants have been experimentally proven to lead to abnormal function. Employing a phenotypic rescue assay with synthetic human messenger RNA variant constructs in shha-/- knockout zebrafish, we evaluated 104 clinically reported in-frame and missense SHH variants. Our data helped us to classify them into loss of function variants (31), hypomorphic variants (33), and nonpathogenic variants (40). We discuss the strengths and weaknesses of currently accepted predictors of variant deleteriousness and the American College of Medical Genetics and Genomics guidelines for variant interpretation in the context of this functional model; furthermore, we demonstrate the robustness of model systems such as zebrafish as a rapid method to resolve variants of uncertain significance.

Rare hypomorphic human variation in the heptahelical domain of SMO contributes to holoprosencephaly phenotypes.

Holoprosencephaly (HPE) is the most common congenital anomaly affecting the forebrain and face in humans and occurs as frequently as 1:250 conceptions or 1:10,000 livebirths. Sonic Hedgehog signaling molecule is one of the best characterized HPE genes that plays crucial roles in numerous developmental processes including midline neural patterning and craniofacial development. The Frizzled class G-protein coupled receptor Smoothened (SMO), whose signaling activity is tightly regulated, is the sole obligate transducer of Hedgehog-related signals. However, except for previous reports of somatic oncogenic driver mutations in human cancers (or mosaic tumors in rare syndromes), any potential disease-related role of SMO genetic variation in humans is largely unknown. To our knowledge, ours is the first report of a human hypomorphic variant revealed by functional testing of seven distinct nonsynonymous SMO variants derived from HPE molecular and clinical data. Here we describe several zebrafish bioassays developed and guided by a systems biology analysis. This analysis strategy, and detection of hypomorphic variation in human SMO, demonstrates the necessity of integrating the genomic variant findings in HPE probands with other components of the Hedgehog gene regulatory network in overall medical interpretations.

BE4max and AncBE4max Are Efficient in Germline Conversion of C:G to T:A Base Pairs in Zebrafish.

The ease of use and robustness of genome editing by CRISPR/Cas9 has led to successful use of gene knockout zebrafish for disease modeling. However, it still remains a challenge to precisely edit the zebrafish genome to create single-nucleotide substitutions, which account for ~60% of human disease-causing mutations. Recently developed base editing nucleases provide an excellent alternate to CRISPR/Cas9-mediated homology dependent repair for generation of zebrafish with point mutations. A new set of cytosine base editors, termed BE4max and AncBE4max, demonstrated improved base editing efficiency in mammalian cells but have not been evaluated in zebrafish. Therefore, we undertook this study to evaluate their efficiency in converting C:G to T:A base pairs in zebrafish by somatic and germline analysis using highly active sgRNAs to twist and ntl genes. Our data demonstrated that these improved BE4max set of plasmids provide desired base substitutions at similar efficiency and without any indels compared to the previously reported BE3 and Target-AID plasmids in zebrafish. Our data also showed that AncBE4max produces fewer incorrect and bystander edits, suggesting that it can be further improved by codon optimization of its components for use in zebrafish.